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Abstract. We present emission line mapping of the strongly barred galaxy NGC 1530 obtained using Fabry-Perot 
interferometry in Ha, at significantly enhanced angular resolution compared with previously published studies. 
The main point of the work is to examine in detail the non-circular components of the velocity held of the gas, 
presumably induced by the strongly non-axisymmetric gravitational potential of the bar. To do this we first 
derive a model rotation curve making minimum assumptions about kinematic symmetry, and go on to measure 
the non-circular component of the full radial velocity field. This clearly reveals the streaming motions associated 
with the spiral density wave producing the arms, and the quasi-elliptical motions with speeds of order 100 km s _1 
aligned with the bar. It also shows in some detail how these flows swing in towards and around the nucleus as they 
cross a circumnuclear resonance, from the dominant "xi orbits" outside the resonance to U X2 orbits" within it. 
Comparing cross-sections of this residual velocity map along and across the bar with the surface brightness map 
in Ha indicates a systematic offset between regions of high non-circular velocity and massive star formation. To 
investigate further we produce maps of velocity gradient along and across the bar. These illustrate very nicely the 
shear compression of the gas, revealed by the location of the dust lanes along loci of maximum velocity gradient 
perpendicular to the bar. They also show clearly how shear, seen in our data as velocity gradient perpendicular 
to the flow, acts to inhibit massive star formation, whereas shocks, seen as strong velocity gradients along the 
flow vector, act to enhance it. Although the inhibiting effect of gas shear flow on star formation has long been 
predicted, this is the clearest observational illustration so far of the effect, thanks to the strong shock-induced 
counterflow system in the bar. It is also the clearest evidence that dust picks out shock-induced inflow along bars. 
These observations should be of considerable interest to those modelling massive star formation in general. 

Key words. Galaxies: general - Galaxies: individual (NGC 1530) - Galaxies: ISM - Galaxies: spiral- Galaxies: 
kinematics and dynamics - ISM: kinematics and dynamics 



1. Introduction 

The interactive relation between dynamics and morpholo- 
gy is a well known aspect of barred galaxies. It has been 
invoked to throw light on one of the more interesting ques- 
tions in galactic evolution: by what mechanism or mecha- 
nisms is material driven in to regions close to the nucleus 
where it can participate in star formation, or contribute 
to the mass of the central supermassive black hole, and 
to the activity around it (Shlosman et al. 1989). As an 
archetypal barred galaxy NGC 1530 has received conside- 
rable attention both morphologically and kinematically. 
These studies have been facilitated observationally by the 
strength and length of the main bar, which compensate for 
its distance. The latter is 36.6 Mpc (Tully & Fisher 1988), 
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which agrees with a more recent calculation of 37 Mpc 
using its recession velocity of 2666 kms _1 (Young et al. 
1989) and a value for H D of 72 kms -1 Mpc -1 (Freedman 
et al. 2001). In spite of this, the bar and its accompanying 
kinematics can be well observed because its deprojected 
radial length is some 14 kpc. The strength of the bar (it 
has one of the highest bar torques measured by Buta et 
al. (2003)) has attracted the attention of observers to the 
distinctive phenomenology of this galaxy. In recent years 
two groups have made major contributions to the study 
of NGC 1530 which are particularly relevant to the work 
presented here. They are Regan et al. (1995, 1996, 1997) 
and the group of Downes & Reynaud (Downes et al. 1996; 
Reynaud & Downes 1997, 1998, 1999; Greve et al. 1999). 
The types of observations these groups present overlap to 
some degree. The group of Downes & Reynaud base most 
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of their work on velocity and surface luminosity mapping 
in CO with some additional stellar information presented 
in Greve et al. (1999), while Regan et al.'s most important 
contributions come from their velocity and surface bright- 
ness mapping in Ha with a useful contribution from CO 
supplemented by near-IR mapping in the initial paper of 
their published series. 

NGC 1530 is a dramatic twin- armed spiral galaxy, with 
arms breaking away sharply from the ends of the long 
primary bar. It is notable that these arms, though trail- 
ing, have shorter opposing lengths, curving inward from 
the tips of the bar as prolongations of the main arms. 
They are well brought out in Ha as seen in Regan et al. 
(1995) or in Fig. of the present paper. The bar ima- 
ged in broad band is wide and contains prominent dust 
lanes which, as the publications by Regan et al. and by 
Downes & Reynaud listed above well point out, pick out 
shocked gas following trajectories generally aligned with 
the bar. These are almost straight along the bar but begin 
to curve inwards from a galactocentric distance of some 
2 kpc, bending sharply in towards the nucleus. They link 
the prominent star-forming regions at the ends of the bar, 
seen very well here in Ha in Fig. with the even more 
prominent star-forming ring feature with radius ^2 kpc 
noted in blue images by Buta (1984) and by Wray (1988) 
and also clearly picked out in Ha in Fig. [fl^ Combining 
near-IR with optical data to give colour maps, Regan et 
al. (1995) explored the stellar and dust morphology in the 
bar and in the central zone following the lanes right into 
the nucleus. One of the aims of the present paper is to 
elucidate further the relation between gas flows and star 
formation which are implied in these morphological data 
described above. 

However, this aim has been pursued with conside- 
rable achievement previously by the two groups men- 
tioned above. Regan et al. (1996) combined atomic hy- 
drogen 21 cm measurements with ionized hydrogen Ha 
Fabry-Perot spectroscopy to derive a rotation curve for 
NGC 1530. From this they derived a bar pattern speed of 
20 kms -1 kpc -1 and used it to predict resonances at radii 
coinciding with prominent morphological features, inclu- 
ding an outer gaseous ring seen in H I, an inner ring seen 
in Hi and Ha, and the nuclear star-forming ring men- 
tioned above. Downes et al. (1996) mapped the galaxy in 
CO and identified the major concentrations of potentially 
star forming molecular gas; these arc located over the cen- 
tral prominent inner disc zone of radius ~2 kpc and extend 
out more weakly over the bar to less prominent concentra- 
tions over the ends of the bar. They added interferometric 
mapping of the central zone, finding a circumnuclear ring 
structure with gas in orbits apparently perpendicular to 
the elongated orbits in the outer parts of the bar. They re- 
late these patterns to the predictions of resonance models 
(Athanassoula 1984; Jenkins & Binney 1994) in which gas 
within corotation moves along X\ orbits along the bar un- 
til it reaches an inner Lindblad resonance (ILR) where it 
moves into the xi orbits perpendicular to the bar. A neat 
schematic diagram in fig. 10 of Downes et al. (1996) illus- 



trates the predicted behaviour of the gas under these cir- 
cumstances. Although the stars follow quasi-elliptical x\ 
orbits, in the case of gas, the shocks are produced along the 
bar and specially at the ends of the bar, where orbit crowd- 
ing occurs, as was well predicted in Roberts et al. (1979). 
However, as the gas reaches the outer limit of the zone 
dominated by X2 dynamics, it suffers another shock, and 
is pulled through a right-angle, falling in again towards 
the nucleus. Gas on X\ orbits with an impact parameter a 
little larger with respect to the nucleus is less affected by 
the effects of the resonance, and proceeds along the out- 
going edge of the bar, albeit with declining velocities. It is 
just this kind of arrangement for allowing gas to fall to the 
centres of barred galaxies which was invoked by Shlosman 
et al. (1989) in their "bars within bars" scenario as a fu- 
eling mechanism for the different types of circumnuclear 
activity: starbursts and AGN. 

Using Fabry-Perot Ha mapping of NGC 1530, Regan 
et al. (1997) demonstrated that the two-dimensional Ha 
velocity field in the bar region shows the skew pattern 
predicted for x\ orbits superposed on the general galac- 
tic rotation. They went on to combine hydrodynamical 
models with their observations to compute a net inflow 
rate to the nuclear zone of ~ IMq per year, with the re- 
maining gas spraying back outwards along the major bar 
axis before coupling into a new inflow cycle. Using molecu- 
lar emission from CO and HCN, the former at angular res- 
olution of 1.8" and the latter at 3.6", Reynaud & Downes 
(1997) detected steep velocity gradients spatially confined 
to two arcs coupling into a ring. These arcs represent shock 
fronts where the inflow along the bar is being constrained 
to change direction and assume the X2 orbital configura- 
tion along the nucleus. These arcs correspond to two small 
arm features picked out in Ha to the north and south of 
the circumnuclear zone of enhanced star formation (see 
Fig. ffli of the present paper). The points where the arcs 
couple to the circumnuclear ring are the main emitters in 
HCN, suggesting points of maximum molecular density. 
Reynaud & Downes postulate the presence of two ILR's: 
one at 1.2 kpc from the nucleus where the shock arcs are 
rooted, and one very close to the nucleus, at the inner edge 
of the molecular ring. They support these inferences using 
resonance theory and their rotation curve, which is de- 
rived kinematically but supported by a dynamical model. 
In a subsequent article Reynaud & Downes (1999) map the 
same area in two emission lines each of 12 CO and 13 CO, 
showing that the molecular gas density is highest in the 
circumnuclear ring and not so high in the arcs referred to 
above. 

Perhaps of most interest and relevance for the study we 
present here is that of Reynaud & Downes (1998), who use 
their CO kinematics of the whole bar region in NGC 1530 
to explore the relation between velocity changes and star 
formation. They conclude that strong velocity changes, 
which represent shocks and shear, are negatively corre- 
lated with star formation and that the distribution of star 
formation does not correspond fully to the distribution 
of available molecular g traced by CO. Although it 
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is true that the major concentration is found in the cen- 
tral circumnuclear zone where star formation is most en- 
hanced, there is considerable CO measured along the bar 
as well as at its ends, while the star formation is clearly 
weak along of the bar. Reynaud & Downes find veloci- 
ty jumps across the dust lanes from velocities of order 
100 kms -1 along the bar in one direction, to similar velo- 
cities in the reverse direction. These jumps maximize half 
way out along the bar from the nucleus on both sides. 
They conclude that steep velocity gradients, both shear 
and shock, prevent molecular clouds from condensing to 
form stars. One of the main aims of the present paper is 
to examine this hypothesis in a little more detail using an 
Ha velocity field. 

We present new high resolution velocity and inten- 
sity maps of NGC 1530, obtained in Ha emission using 
a Fabry-Perot interferometer. The instrument provides a 
three-dimensional "data cube" , akin to a conventional ra- 
diastronomical set of line maps, in which the emission in- 
tensity is given as a function of three coordinates: two 
positional coordinates and the velocity. From this "data 
cube" one goes on to compute maps in integrated line 
intensity, line peak velocity, and velocity dispersion. The 
next step is to obtain the rotation curve of the galaxy 
and to subtract it off two-dimensionally from the com- 
plete observed velocity field. This allows us to obtain the 
non-circular motions due to non-axisymmetric elements of 
the galaxy. For NGC 1530 the dominant source of non- 
axisymmetry is the strong bar of length almost 28 kpc, 
which dominates the appearance of the galactic disk. We 
also find kinematic signatures of the streaming motions in 
the sharply defined spiral arms. Besides we also see, in the 
central 5 kpc, evidence for a distinct kinematic structure 
predicted by classical resonance theory, and already care- 
fully observed and characterised in the studies by Regan 
and coworkers, and by Reynaud & Downes in the referen- 
ces given above. 

As a further step we have prepared plots of the velocity 
gradients in the whole of the bar and arm structure as ima- 
ged in Ha emission. We have related the velocity gradient 
morphology to the distribution of the massive star forma- 
tion in arms, bar, and the central zone. Our main steps 
forward compared with previous studies are due to the 
uniformly high angular resolution we could obtain over the 
full image. This has enabled us to make a more precise de- 
finition of the rotation curve, and hence an improved map 
of the non-circular velocity field, both around the main 
bar, and in the central zone. From the two-dimensional 
display of the non-circular velocities we have gone on to 
derive a map of the velocity gradients, with a resolution 
better than 2" over the whole field of the galaxy. Such a 
map has not been previously produced, so we are able to 
make progress here in characterising the geometrical and 
physical relationships between velocity gradients and star 
formation. 



2. Observations and data reduction 

The Fabry-Perot observations of NGC 1530 were taken 
at the Roque de los Muchachos Observatory with the 
4.2m William Herschel Telescope (WHT) on the night 
of September 2nd, 1998. The TAURUS-II instrument, 
equipped with a TEK CCD camera was used to observe 
the galaxy in the Ha emission line, with a plate scale 
of 0.292" per pixel, which corresponds to a pixel scale 
179 pc arcsec -1 at the distance of the galaxy (37 Mpc). 
The observations consisted of exposures of 150 s at equally 
spaced positions of the etalon, at intervals of 0.407 A, co- 
rresponding to 18.62 kms -1 . The total number of steps 
was 55, allowing us to scan the full wavelength range of 
the Ha line in the galaxy. We used the appropriately red- 
shifted narrow band Ha filter as an order-sorting filter 
(A C ~6613A, AA=15A, corresponding to the galaxy's sys- 
temic velocity v sys =2461 km s -1 ; de Vaucouleurs et al. 
1991; hereinafter RC3). Wavelength and phase calibration 
were performed using observations of a calibration lamp 
(taken before and after the observations of the galaxy) , u- 
sing the TAUCAL package in the FIGARO environment. 
The spatial resolution of the resulting data set is The 
data set was placed on a correctly oriented spatial grid 
by comparing positions of point-like features in the Ha 
Fabry-Perot intensity map with those in the narrow-band 
Ha image from Relaho et. al (2003b), yielding maximum 
uncertainties in positions of 0.3". The final calibrated data 
cube has 941x941 pixels x 55 'planes' in wavelength, se- 
parated by 18.62 kms -1 . We examined plane by plane the 
data cube and found that it had a considerable number 
of cosmic ray events. We eliminated them with the task 
imedit in IRAF, which allowed us to replace the anoma- 
lous pixels by a mean background value. We determined 
which channels of the data set were free of Ha line emission 
after smoothing the data to a resolution of 2.5" x 2.5". We 
found 10 channels free of line emission on the low- velocity 
side, and 13 channels on the high-velocity side. The con- 
tinuum was determined by fitting a linear relation to the 
line-free channels, and was then subtracted from the line 
emission channels. 

In order to study the kinematics on different scales, the 
data-cube was convolved down from the original resolu- 
tion of l"xl" to resolutions of 2"x2", 3"x3", 6"x6"and 
10"xl0". For each data set we obtained the total Ha 
intensity (zeroth moment), velocity (first moment) and 
velocity dispersion (second moment) maps. The proce- 
dure to obtain the moment maps is equivalent for each 
smoothed data cube so we describe the procedure followed 
for the full resolution data set only. A reader wanting a 
fuller description of the procedure can find it in Knapen 
(1998). The points with low signal were eliminated pro- 
ducing a conditionally transferred data cube, in which 
values were retained only at positions where the inten- 
sity in the smoothed data cube of 3" x 3" was larger than 
2.5 times the rms noise of the smoothed cube; pixel va- 
lues at all other positions were set to undefined. We then 
removed noise peaks outside the area where Ha emission 
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Fig. 1. a) Intensity map (zeroth order moment) of the Ha emission in NGC 1530. The map was obtained from the 
2" resolution data cube. The kinematic centre is marked with a black cross, b) Velocity map (first moment) of the 
Ha emission in NGC 1530, derived from the 2" resolution data cube, c) Model velocity map of NGC 1530, obtained 
from the Ha rotation curve of the galaxy shown in Fig. [21 with Ha isointensity contours overlaid. We can clearly 
see "wiggles" at the radii of the spiral arms, due to the density wave streaming motions not subtracted off from 
the rotation curve, d) Map of the residual velocities of NGC 1530, obtained by subtracting the model velocity field 
(Fig. ffl;, after taking out the ripples due to the streaming motions in the spiral arms, see text for details) from the 
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is expected by setting the pixel values there to undefined 
too. This was done interactively by inspecting the high 
resolution channel maps one by one and comparing them 
with the corresponding and adjacent channel maps in the 
3"x3" data set. The resulting data cube was then used to 
calculate the moment maps using the moments task in 
the GIPSY package. 

moments obtains the zeroth, first and second mo- 
ments of the emission at each spatial pixel in the galaxy 
and produces a corresponding map for each moment. In 
the procedure, moments takes into account the signal 
only when two conditions are met: the signal is higher 
than 2.5 times the rms noise of the high-resolution cube, 
in at least three adjacent channels. The intensity map for 
the 2" resolution data cube is shown in Fig. and the 
first moment (velocity) map is shown in Fig. ^d. 

Although the moments procedure is a more useful tech- 
nique for studying the global kinematics in the galaxy, it 
must be treated with caution when it is used for a detailed 
description of the kinematics of individual Hn regions. 
The moments procedure gives an underestimate of the ve- 
locity dispersion of the profile (van der Kruit & Shostak 
1982), thus, if a precision measurement of this quantity is 
required, a better fit to the line profiles, using Gaussian 
functions, is needed. 
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Fig. 2. Ha rotation curve of the disc of NGC 1530, for 
the receding and the approaching sides of the galaxy 
separately. The agreement between the two sides shows 
the symmetry of the velocity field and the reliability of 
the rotation centre. 



3. Kinematics 

3.1. Rotation curve 

The rotation curve of the galaxy can be derived from the 
velocity map, which gives the projected velocity along the 
line of sight in each position of the galaxy. Assuming that 
the measurements refer to a position on a single inclined 
plane, the projected velocity can be related to the rotation 
velocity in the plane of the galaxy by this expression: 

v(x,y) = v sys + « rot (r)cos(0)sin(i) (1) 

where v(x,y) denotes the radial velocity at the sky coor- 
dinates x and y, v sys the systemic velocity, w rot the ro- 
tational velocity, i the inclination angle of the normal to 
the galaxy plane to the line of sight, and 9 the azimuthal 
distance from the major axis in the plane of the galaxy, re- 
lated to the rotation centre (xq, yo) and the position angle 
of the galaxy, PA. We derived the rotation curve following 
the procedure described by Begeman (1989), where the 
rotation curve is determined by fitting tilted rings to the 
observed velocity field. For each ring independently, the 
parameters w sys , V I0 t, % and Xq, yo an d PA are obtained 
making a least-squares fit to Eq. using the observed 
velocity field v(x,y). 

We discarded the points within a certain free-angle 
$max around the minor axis since they carry less informa- 
tion about the underlying circular velocity of the galaxy 
and the values obtained in this zone are less accurate due 
to projection effects. In the case of NGC 1530 a value of 
#max = 45° was used to avoid the non-circular motions 
along the bar. Besides, we give more weight to positions 



close to the major axis using a weighting function pro- 
portional to | cos 8\ in the fit. The procedure used to ob- 
tain each of the parameters which characterise the rotation 
curve is the following. 

The parameters v sys , xq and yo depend on the symme- 
try of the velocity field and are the same for each ring. 
We obtained the position of the kinematic centre (xn, yo) 
using the velocity map of 6" resolution. We took annuli of 
3" width, fixed values for the heliocentric systemic velo- 
city at 2461 km s -1 (RC3), inclination angle at i = 50° (a 
mean between the values of Downes et al. 1996 and Regan 
et al. 1996) and position angle at 188° (the value found by 
Regan et al. 1996) and obtained (xo,yo) using Eq. ^ We 
found a large scatter in the values of the kinematic centre 
for annuli far from the centre of the galaxy. This is due to 
the paucity of Ha emission in the disc of the galaxy outside 
the arms and to the fact that lower resolution data cubes 
are more spatially affected by the smoothed non-circular 
motions along the bar in the data cubes. The mean value 
of (xq, yo) for the rings inside 12" corresponds to the coor- 
dinates 4 h 23 m 26.94 s (±0.02 s ) 75° 17' 44.0"(±0.1"), J2000. 
We checked this value using the same procedure for the 
velocity maps of lower and higher (10" and 2") resolution, 
and found that the differences in the rotation centre from 
the value adopted were inside the error bars. 

The systemic velocity was obtained from the 10" reso- 
lution velocity map with rings of width 5". The centre po- 
sition was fixed to the value found previously, and the in- 
clination and position angle were fixed to the same values 
used to derive the rotation centre. The mean value for v sys 
in all the rings is 2466 ± 4 kms -1 , where the error is the 
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standard deviation from the mean. This value is in agree- 
ment with the value found by Reynaud et al. (1997) of 
u sys = 2466 kins -1 for CO data, but differs from that ob- 
tained by Regan et al. 1996 in Ha, v sys = 2447±2 kms -1 . 
A way to test whether the rotation centre and the sys- 
temic velocity are well determined is shown in Fig. [3 
where we have plotted rotation curves for the receding 
and approaching halves of the galaxy derived separately u- 
sing the rotation centre and systemic velocity found above. 
The coincidence of the velocity points for both sides, espe- 
cially in the central zone of the galaxy, provides evidence 
of the symmetry of the velocity field and the validity of 
the values for the rotation centre and systemic velocity. 
We also derived the rotation curves for the receding and a- 
pproaching halves with the rotation centre found in Regan 
et al. (1996), but in that case we found discrepancies which 
appeared unphysically large. It is the better resolution of 
our data (1", while Regan et al. (1996) has a net resolution 
of 4") which allows us to obtain a more accurate position 
for the rotation centre. 

Finally, we obtained the rotation curve using the ve- 
locity map at 2" resolution in Fig^>. For the internal part 
of the galaxy, out to 12.5", we used rings of width 1", and 
rings of width 6" for larger radii. The rotation curve ob- 
tained in this way is shown in Fig.|3 where the PA and the 
inclination were fixed at 188° and 50° respectively. Fig. 
also displays the rotation curve from Regan et al. 1996, 
which is similar to that found here, but we obtain higher 
values of the rotation velocity in the inner part, 10"-15", 
and also in the 70"-90"radial range. In order to check the 
PA adopted, we left it as a free parameter in the fit, find- 
ing that the values given by the program for most of the 
annuli are close to 188° and that the rotation curve did not 
change significantly from the curve with a fixed PA. The 
adopted inclination angle i = 50° was also checked. We 
fitted w ro t and i with Eq. ^ keeping the other kinematic 
parameters fixed, and found a mean value of i — (53 ±8)°, 
in good agreement with the value adopted. 

To test the importance of a careful definition of the 
the free-angle we reduced its value to 20° and rederived 
the rotation curve. As seen in Fig. 0] this produced "ro- 
tational velocities" rising to values of over 500 kins" 1 , as 
there is an obvious contribution from the non-rotational 
motion along the bar projected into the line of sight. As we 
will see below the amplitude of the gas motion along the 
bar is over 100 kms -1 , and it is easy to see that this will 
yield the unrealistically high velocities shown in Fig. 0] 
Fortunately, the bar orientation in NGC 1530 does not 
make it too difficult to isolate and avoid the emission from 
gas which contains a significant non-circular velocity com- 
ponent, but Fig. goffers a warning to those using maps of 
much lower resolution in similar circumstances. 



3.2. The residual velocity map 

After deriving the mean rotation curve, an easy way to 
explore non-circular motions within the galactic disc is 
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Fig. 3. Ha rotation curve of the disc of NGC 1530, up 
to 18 kpc from the galactic centre. The rotation curve 
was obtained by fitting tilted rings to the observed ve- 
locity field in Fig. ^J). Also shown is the rotation curve 
obtained by Regan et al. (1996). 



to make a two-dimensional projection of this curve and 
subtract it from the observed velocity field. The two- 
dimensional projection of the initial rotation curve is seen 
in Fig.^;, where we have superposed an isophote map of 
intensities to show the positions of the emitting compo- 
nents. The projection has classical symmetry, but shows 
a few small features: "wiggles" which occur at some three 
quarters of the radial distance from the centre to the edge 
of the field. These can be explained as effects of stream- 
ing motions in the spiral arms, which have not been sub- 
tracted adequately from the mean rotation curve. 

While it is difficult to get an exact separation of ro- 
tation and streaming motions, we have made a first a- 
pproximation by assuming that the true rotation curve 
is smooth. Applying a second order polynomial fit to the 
part that is observed as smooth, i.e. between 27" and 63" 
from the centre, we then extrapolate this to cover the 
range between 63" and the edge of the field at 99" from 
the centre, and go on to combine it with the original curve 
out to 63" radius, to give the final "true" rotation curve. 
The result of this procedure is shown in Fig. where we 
have compared our original rotation curve with the new 
curve and also with the curve by Regan et al. (1996). Our 
procedure has isolated the wiggles between 63" and 93". 

We can then proceed to produce a new two- 
dimensional projected model for the rotation in which the 
wiggles are no longer present, and to subtract off this new 
rotation model from the original velocity field, yielding a 
resulting map of residuals as shown in Fig. QJ1. There are 
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Fig. 4. Ha rotation curves of the disc of NGC 1530 ob- 
tained with different "free-angles" and with rings of 1" 
width. For free-angles less than ~ 30° there is an obvi- 
ous contribution from the non-rotational motion along 
the bar in the mean calculated rotational velocity that 
can "contaminate" the rotation curve of the galaxy. 



Fig. 5. Ha rotation curve of the disc of NGC 1530. This 
rotation curve has been obtained by subtracting the ef- 
fect of the streaming motions in the spiral arms of the 
galaxy from the mean rotation curve shown in Fig. 
(see text for further details). 



clearly ordered fields of non-circular velocity both around 
the bar, and in the arms, but we can check on the mo- 
del subtraction by examining the residuals of the regions 
which are in the disc outside the arms and bar. These in 
general show, as they should if the rotation curve sub- 
traction is valid, net residual velocities close to zero: less 
than ~10 km s -1 , which is close to the precision limit of 
the observations. We can now see clearly the streaming 
motions across the arms in Fig. ^L. These have projected 
values between 10 and 30 km s _1 , with similar amplitudes 
in both arms. 

3.2.1. The residual velocities along the bar 

Of most striking interest is the fairly complex non-circular 
velocity field associated with the bar. Approaching veloc- 
ities of up to 120 km s _1 are observed along the north 
side of the bar, and receding velocities of the same am- 
plitude along the south side. These values, and the im- 
plied vectorial direction in elongated quasi-elliptical or- 
bits around the bar with long axis parallel to the bar axis, 
are in line with classical predictions of the dynamical ef- 
fects of a stellar bar on the rotating gas which date back 
to Huntley et al. (1978) and Roberts et al. (1979) as well 
as to more recent studies by Athanassoula (1992). Most of 
the observational work has been performed either in H I or 
more recently in CO emission. There have been relatively 
few Ha studies of this phenomenon, for two probable rea- 
sons: slit spectroscopy is costly in time if more than a 
couple of position angles are taken (one good early study 



in this mode using six slit positions on NGC 1512 was 
published by Lindblad & Jorsater 1981); and Ha emis- 
sion is much patchier than either Hi or CO. However, 
as we can see from the data in the present article, the 
use of a scanning Fabry-Perot for this work permits full 
two-dimensional coverage at an angular resolution better 
than that of essentially all H I or CO measurements, which 
can compensate to an interesting degree for the intrinsic 
patchiness of the observed field. 

In Fig. there are several complexities in the non- 
circular velocity field which are worth pointing out. Taking 
the north side of the bar we can see a zone with high bar- 
streaming velocities, of order 120 km s _1 in a direction 
approaching us to the east of the centre, and another re- 
gion, more restricted in size and velocity range, with simi- 
lar velocities to the west. To the south of the bar we can 
see a diametrically symmetric flow, with high receding ve- 
locities, again of order 120 km s _1 receding on west side, 
and two smaller zones with similar velocities to the east. 
It is tempting to consider these flows as direct evidence for 
motion aligned with the bar and terminating in zones of 
shock towards the ends of the bar. We will see later that 
there is considerable evidence in favour of this, and of the 
direct causal links between properties of these flows, and 
both the stimulation and inhibition of massive star forma- 
tion. 

Examining Fig. QJi, the question arises as to how the 
flow can apparently be stopped halfway along the bar, yet 
continue on beyond the central zones as described above. 
This appearance does not give a true picture of the ve- 
locity field and is essentially due to the patchiness of the 
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emission. However taking the latter into account we can 
interpret the observations in terms of a resonance model 
of the type illustrated in Downes et al. (1996) to interpret 
their measurements in CO, based on models incorporating 
an ILR to a main galactic bar by Athanassoula (1992). 
Qualitatively, these models show how the bar dynamics 
produces orbits termed Xi, whose major axes are aligned 
with the major axis of the bar, but that the presence of 
an inner Lindblad resonance within the bar structure (not 
far from the nucleus) implies that gas will move on the 
perpendicular X2 elliptical orbits within the ILR. The be- 
haviour of a particular gas cloud in an x\ orbit will then 
depend critically on its impact parameter with respect to 
the nucleus, i.e. on the length of its orbital semi-minor 
axis. Gas on orbits close to the bar major axis will be 
dragged into X2 orbits as it passes close to the nucleus, 
and will not escape along the bar. However, gas on orbits 
well away from the major axis will remain essentially on 
those orbits until it reaches the far end of the bar, where 
it will be braked and shocked. Gas on intermediate orbits 
will swing in towards the bar major axis, where in turn it 
will also be braked and shocked. In Fig.EJwe have illustra- 
ted these types of orbital behaviour, with a representative 
set of five trajectories on either side of the bar. Gas mo- 
ving on the innermost trajectory will be pulled into an x-i 
orbit, while gas on the others will end up impinging on the 
bar progressively further along, as their initial impact pa- 
rameters increase. Because the velocity vector is observed 
only where star formation is occurring, we cannot follow 
these trajectories along their lengths, but see clearly their 
initial and final points in the residual velocity field. 

3.2.2. Residual velocity cross-sections along and 
perpendicular to the bar 

In Figs EJi, Et> , Et and EH we show profiles in residual 
velocity perpendicular and parallel to the bar direction 
(PA=116°). The positions and directions of the cross- 
sections chosen are plotted over the residual velocity map 
in Fig For profiles perpendicular to the bar and loca- 
ted where the bar connects with the spiral arms (-40" in 
Fig. Et and 40" in Fig. Et>) the absolute value of the re- 
sidual velocity along the line of sight is ~70-80kms _1 . 
The velocities rise as we approach the centre of the bar 
and at ±30" are ~100kms _1 . We note particularly the 
symmetry of the profiles on cither side, at ±20" of the 
profile through the nucleus. We can also see very clearly 
the steep velocity gradients at these positions. In the pro- 
files parallel to the bar (FigEt and FigEtO there is also 
evidence of rapid velocity changes. In particular the pro- 
files at 4" and 8" in Fig El:, and -12" and -8" in Fig EH, 
show very sharp velocity gradients. These gradients mark 
a combination of braking and projection effects: braking 
effects as the gas flow along the bar meets the central 
gas concentration, and projection effects as the flow vec- 
tor switches direction from nearly along the line of sight 
(along the bar) to nearly perpendicular to the line of sight 



(as the gas is pulled into circumnuclear orbits). A result of 
examining these profiles was to induce us to make a more 
detailed map of the velocity gradients, which we will show 
and discuss below. 

3.3. Residual velocity gradient maps 

Following up the indication of the presence of steep ve- 
locity gradients in Figs. Et, Et>j Et an d Eli, we opted to 
produce a full map of the velocity gradients in the sys- 
tem, taking advantage of our excellent angular resolution. 
To make sure that the steep gradients were not artefacts 
of our data processing we first used the 2" resolution data 
cube, and made position-velocity diagrams at those points 
where Figs. Et,Et>,Et, Eli indicated that we ought to find 
rapid velocity changes. This exercise was non-trivial, since 
the steepest gradients occur in zones of low Ha intensity, 
for reasons which we will discuss further below. The dia- 
grams are shown in Fig. 1101 where we can see at once that 
the velocity gradients are real and very steep, at the posi- 
tions chosen. We go on to derive maps in velocity gradient, 
selecting two directions, parallel and perpendicular to the 
bar, to optimize the information about the gas flows in- 
duced by the gravitational perturbation of the bar. The 
gradient for a given direction was derived in a straight- 
forward way as follows: the velocity residual map was 
displaced by one pixel in that direction, and subtracted 
off from the original map; this procedure was repeated 
in the opposite direction, and the gradient was found by 
co-adding the absolute values for the two maps, and nor- 
malizing the result, dividing by two. The same procedure 
was followed for both of the selected directions, yielding 
two independent velocity gradient maps. The whole map 
of residual velocity gradient in the direction perpendicular 
to the bar is shown in Fig.Et, while in Figs. II Hi and II lb 
we show maps of the gradients in residual velocity, mea- 
sured in the directions perpendicular to the bar and along 
the bar respectively, in the zone of the bar. The regions of 
steep gradient show up very clearly, and reach high val- 
ues (~ 0.35 kms _1 pc~ 1 ). Near the centre of the bar we 
can pick out elongated regions where the motions along 
the bar axis are in counterflow, so that an amplitude of 
over 100 kms -1 in one direction rapidly gives way to the 
same amplitude in the opposite direction. The measured 
gradients imply that this occurs over distances of a few 
hundreds of parsecs, close to the limit of direct detection, 
so that there may in fact be much steeper shocks than 
this. Our decomposition of the gradients into these two 
perpendicular directions does not effectuate a very clean 
separation between different intrinsic dynamical effects, 
above all because the observations give the projected com- 
ponents along our line of sight, and projection effects are 
not removed. This would be difficult to do correctly with- 
out a previous underlying model, which would assume too 
much about the directions of motion to be reliable, so we 
prefer to leave the observations as they are. In general, 
zones of high velocity gradient in one direction also show 
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Fig. 7. Lines tracing the cross-sections perpendicular and 
parallel to the bar of NGC 1530 represented in Figs. 
to|SJi over a portion of the residual velocity map. 




Fig. 8. Schematic representation of the gas orbits along the 
bar on the residual velocity map, illustrating the depen- 
dence of the gas trajectories on orbit semi-minor axis. 
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Fig. 9. a) Velocity gradient map perpendicular to the bar. Zones of maximum shear are clearly seen (cf. the dust 
lanes in Fig. Eh)- b) Image of NGC 1530 in the optical I^band obtained with the KPNO 4 meter Mayall telescope 
(NOAO / AURA/NSF) . 
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In high gradients in the perpendicular direction, though 
the details differ on very small scales. Fig. 03 we have 
also included an image of NGC 1530 in the Fband. The 
morphological and positional agreement between the dust 
lanes and the zones of maximum velocity gradient perpen- 
dicular to the bar illustrated in this figure is striking. This 
result may not be qualitatively unexpected, since various 
models of gas flows and shocks in barred galaxies do pre- 
dict that gas and dust densities should be elevated in the 
shocked zones. However, to our knowledge, this is the first 
presentation of a direct one-to-one mapping of the veloc- 
ity gradients in shocks and the dust lanes along a bar with 
comparable angular resolution. 

Comparing these velocity gradient maps, and in par- 
ticular the map of gradients perpendicular to the bar, 
Fig-EK orlllH. with the total intensity map in Fig. we 
can see clearly that the star formation is peaking in zones 
where the velocity gradients are low. There is a very good 
anticorrelation here, which holds over the whole of the face 
of the galaxy: in the arms, along the bar, and in the central 
zone. In general, when comparing the non-circular veloci- 
ties in Fig. Hlb with the velocity gradients in Fig.lllb.d. we 
find high velocity gradients associated spatially with sig- 
nificant non-circular velocities, but there are two specific 
zones where this is clearly not the case. Along the NW side 
of the bar there are two bright star-forming regions, the 
less intense of which is situated within gas with a projected 
residual flow vector of amplitude ^100 kms" 1 , as seen in 
Fig. II lb . This same zone, is clearly in a velocity gradient 
minimum as seen in Figs. 1111 1 and lllb . A second exam- 
ple is the intense star forming zone where the bar meets 



the NW arm, where the velocity gradients are low, but 
the net velocity vector has an amplitude of ~30 kms -1 . 
Both of these are regions with strong star formation, and 
low velocity gradient but relatively high residual velocity. 
We will discuss the implications of these correlations and 
anticorrelations in Sect. 4. 

3.4. The velocity dispersion 

The second moment of the intensity-wavelength distri- 
bution represents the velocity dispersion of the emission 
line observed. Since we are interested in the dispersion 
due to the turbulent (non-thermal) motions of gas, we 
must extract it from the observed Ha emission line pro- 
files taking into account effects which also contribute to 
broaden the observed line shape. These are the natural Ha 
emission line broadening, thermal broadening and instru- 
mental broadening. We subtracted these components in 
quadrature from each point of the observed velocity dis- 
persion map, to obtain a new map of the dispersion due 
to the non-thermal motions of the ionized gas. For further 
details on the procedure see Rozas et al. (1998, 2000). In 
Fig. II lb we show the map of non-thermal velocity disper- 
sion in the zone of the bar. In general there is a close 
correlation between high surface brightness (as seen in 
Fig. ll"Tk ) and high velocity dispersion, as expected know- 
ing that much of the measured turbulence is stimulated 
by the energy output of young massive stars, which also 
give rise to the Ha. There are some zones, notably towards 
the ends of the bar, where this agreement is less clear, i.e. 
where the ratio of the Ha surface brightness to the mea- 
sured dispersion is much lower than for other star forming 
zones. We attribute this to the fact that our measurements 
of dispersion cannot avoid including components of steep 
localized velocity gradients which are particularly impor- 
tant where the arms break away from the bar. As we will 
see in Sect. 14.21 there is a clear anticorrelation between the 
local star formation rate and the gas velocity gradient, and 
the observations of dispersion alone would not easily pick 
this up. 

4. Relationships between the velocity field and 
the local star formation rate 

4.1. Bar cross-sections in velocity and Ha surface 
brightness 

To examine more closely the relation between the current 
star formation rate as measured in Ha surface brightness, 
and the velocity field of the gas in the region of the bar, 
we obtained a series of cross-sections in which the resi- 
dual velocity can be compared directly with the Ha sur- 
face brightness, as displayed in Fig. IT5K. and in Fig. H"2"b 
with the velocity dispersion. In Fig. 112b we show a set 
of these cross-sections perpendicular to and parallel to 
the bar, comparing surface brightness with residual velo- 
cities, i.e. non-circular velocity. As a general feature of the 
cross-sections it is clear that there is an anti-correlation 
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between Ha intensity and residual velocity 1 . We can make 
a plausible hypothesis that suggests two distinct pheno- 
mena, one of which inhibits star formation and the other 
of which helps to produce it, which combine to yield this 
general result. Around the bar, the long elliptical orbits 
should have (see e.g. Huntley et al. 1978) velocity vec- 
tors around the ellipse which are high close to the bar 
major axis, and fall further and further from the axis. In 
this configuration the maximum residual velocity around 
the bar should be close to the zone of maximum veloc- 
ity shear. The conditions which determine the stability of 
clouds for global star formation in galaxy disks have been 
shown by Kennicutt (1989), following Toomre (1964), to 
depend on the relationship between the gas surface den- 



1 To make for a homogeneous comparison we have plotted as 
positive the amplitude of the velocity vector without reference 
to its direction, while from the map in Fig.^i it should be clear 
that the velocities are in opposite directions along opposite 
sides of the bar. 



sity and the velocity dispersion in the plane. However, 
this general condition should be affected by the presence 
of velocity shear, which would act against the condensa- 
tion of massive gas clouds to form stars, as suggested by 
Reynaud & Downes (1998). Therefore, in the global mo- 
tion around the bar we might expect to find a lower star 
formation rate in zones of high velocity. However, within 
this global framework, certain shocked zones are predicted 
to be places where star formation is accelerated (see e.g. 
Elmegreen et al. 2002), and where the global motion is lo- 
cally strongly braked. This will give rise to zones of low ve- 
locity and high star formation rate which also contribute 
to a negative correlation between gas velocity and star 
formation rate, but with variations over more compressed 
spatial scales. The global anti-correlation is well seen in 
all the cross-sections of Fig. 112b . where peaks in the Ha 
surface brightness nearly all coincide with troughs in the 
residual velocity. Local effects are especially well seen in 
the cross-section through the nucleus parallel to the bar, 
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which cuts two shock fronts; near these the star formation 
rate is low, but in the radial galactocentric range between 
them it is very high. These relationships are also seen in 
Fig. IT2b . where we have made cross-sections along the 
same tracks, perpendicular to and parallel to the bar, but 
where we have plotted the velocity dispersion instead of 
the Ha surface brightness. Qualitatively the graphs are 
very similar, with a similar set of anticorrelation patterns. 
The main difference is that the velocity dispersion peaks 
tend to be broader than those in surface brightness, re- 
flecting the turbulent velocity structure within individual 
ultraluminous Hil regions (there is a clear tendency for 
the brighter regions to have the biggest turbulent widths; 
see e.g. Rozas et. al 1998; Relaho et al. 2003b, in prepa- 
ration). We should also note velocity gradients on diffe- 
rent scales, in particular gradients in the global velocity 
pattern can contribute to an observed local value of the 
velocity dispersion, while there is no analogous effect for 
surface brightness. There is quantitative material in these 
data which can contribute to theoretical modelling of sce- 
narios for massive star formation, though this is beyond 
the direct scope of the present article. 



4.2. Bar cross-sections in velocity gradient and Ha 
surface brightness 

We saw qualitatively, by inspection in Sect. 13.31 that where 
the velocity gradients in the bar are high, star formation 
appears less intense. For a more quantitative appreciation 
of this effect we have made graphs of residual velocity gra- 
dient and Ha intensity, taken along scans perpendicular 
to and parallel to the bar. We have used both Figs. II lb 
and II lfc l for this, so that we are comparing Ha surface 
brightness with velocity gradient parallel to the bar and 
perpendicular to the bar, respectively in the two cross- 
sections in Figs.lT3"h andlT3b. 

In both graphs and in all the cross-sections, we note 
the coincidence between the intensity peaks and local 
minima, or points of inflection, in the velocity gradient. 
There is a striking anticorrelation between velocity gra- 
dient and star formation rate as represented by the Ha 
surface brightness, and it is particularly striking in zones 
of high velocity gradient. In the two zones where the ve- 
locity gradient is both high and rapidly varying, between 
10" and 15" from the centre of the scans parallel to the 
bar through the nucleus (centre-right plot of Figs. I13h 
and 113b). there is a notable absence of star formation. 
These are the shocked zones at the edges of the central Xi 
dominated region, where the flow along the bar is being 
braked and the gas is changing its orbital configuration. 
The star formation within the central region is occurring 
well within this radial range, i.e. well offset from these 
shocks. As well as the global effects just described, there 
are local dips in the velocity gradient, of smaller ampli- 
tude, coincident with the most intense star forming re- 
gions, which reflect the kinematic impact of the gas mo- 



tions within the individual regions on the overall velocity 
field. 



5. Velocity distribution and morphology within 
the inner disk 

5.1. The nuclear bulge and the circumnuclear velocity 
pattern 

Both morphologically and kinematically the central zone, 
which shows up as a projected disk of major axis radius 
~ 10", i.e. ~ 1.8 kpc, is easily distinguished in Ha from 
the bar. At its centre, as well seen in the Hubble Space 
Telescope (HST)-mCMOS image presented in Fig. Et, 
there appears to be a very small bulge, less than 2" in 
diameter, detectable via its circular isophotal shape. The 
rotation curve for the galaxy as a whole, which we pre- 
sented in Fig- EI shows an apparently enhanced rotational 
velocity gradient within the inner ~ 5", but we had pro- 
duced this curve on the assumption that the only signifi- 
cant velocity component here is in the tangential direction 
to circular orbits in the plane of the galaxy. 

From Fig. 114b we find that there is no major extra 
mass component on the scale of this steeper gradient, so 
we explored the possibility that it is in fact indicating the 
presence of non-circular motions within this inner disk. To 
do this we assumed, from the form of the rotation curve 
as a whole, that the velocity gradient between 5" and 10" 
represents the true rotation curve, and can be interpo- 
lated between this portion of the curve and the nucleus. 
We then subtracted off the rotated interpolated curve, to 
produce a new residual velocity pattern for the central 
disk, which is shown in Fig. 114b . We can see from this fi- 
gure that as well as the previously identified non-circular 
velocities beyond 10" from the nucleus, we now have two 
lobes of oppositely directed velocity, nearly symmetric, to 
the north and south of the nucleus, with peak projected 
values of ^50 kms -1 at a little under 5" from the cen- 
tre. This pattern might be attributable to an inner bar or 
oval distorsion within this inner disk. To check this out 
we made an unsharp masked image (see Erwin & Sparke 
2002) from the NICMOS near IR image shown in Fig. HIL 
by subtracting off from the original a version smoothed 
with a gaussian of 5 pixels half-width (~ 0.38" FWHM). 
The result is shown in Fig. 1151 in which we can see a 
bright nucleus, and a small asymmetric spiral, possibly e- 
merging from a tiny bar, whose major axis appears to be 
at PA^ 120°, but which is not well enough resolved even 
at this resolution. 

The dynamical scope of this structure can be appre- 
ciated from Fig. 114b . where it has been plotted onto the 
residual velocity map of the inner disk, and where we can 
see that it is far too small to be causing the lobed velocity 
structure, which demands a different explanation. 
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Fig. 12. a) Profiles of the normalized surface brightness in Ha emission and residual velocity parallel and perpendicular 
to the bar. The distance of each profile track from the kinematic centre is indicated in the top right corner of each 
plot. Positive distances in the perpendicular and parallel profiles are respectively towards the NW and NE in Fig. 
b) Profiles of the absolute residual velocity and nonthermal velocity dispersion parallel and perpendicular to the bar. 
The tracks of the profiles are the same as those in Fig. 112b . 



5.2. The inner kinematic structure: a spiralling inflow 
pattern 

We can go on to interpret the inner residual velocity 
structure as follows. The gas flow along the bar and 
its deviation into circumnuclear orbits can be followed 
via the presence of dust lanes which indicate compressed 
and shocked gas. These dust lanes in the central disk of 
NGC 1530 were imaged by Perez-Ramirez et al. (2000; 
their fig. 2b) in the near IR, at subarcsec resolution, and 
are especially well seen in their J-K colour image, which 
we reproduce here in Fig 114b . Two clear dust lanes can 
be seen to the north and south of the inner disk, coming 
inwards along the direction of the bar major axis, then 
curling around in orbits which take them to the far side 
of the nucleus in both cases, where their trajectories have 



swung through 270°, and terminating on opposite sides of 
the nucleus at some three arcsec from the centre, along 
a line with PA~120°. The positions of the centres of our 
flow lobes in Fig. 114b coincide with the inward curving 
dust lanes, so we attribute them to gas flow coinciding 
with these lanes. The observed velocity structure of the 
lobes is strongly affected by the fact that we observe the 
velocities along a single line of sight, and so we detect 
only their projected components. Where the velocity vec- 
tor of the inflowing gas is perpendicular to the line of sight 
the flow cannot be detected. Thus, the lobe to the north 
of the nucleus is in fact reflecting the final stages of the 
inflow from the southern edge of the bar, and vice versa. 
The two flows spiral around one another in an interlocking 
pattern. As a test of this prediction we have produced a 
simple kinematic model showing how these flows would 
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Fig. 13. a) Profiles of the normalized Ha surface brightness and residual velocity gradient parallel to the bar. The tracks 
of the profiles are the same as those in Fig.ll2fa. b) Same as (a) but for the residual velocity gradient perpendicular to 
the bar. Note the offsets and anticorrelations of the intense star forming regions and regions of high velocity gradient 
shown in this figure and in Fig.H3h (see text). 



appear in projected velocity. We have assumed for sim- 
plicity that the intrinsic velocity vector has a constant 
amplitude, equal to 60kms _1 , and projected it into the 
line of sight to the observer, using the inclination and posi- 
tion angle of the galaxy, and assuming that the dust lanes 
in the observations from Perez-Ramirez et al. (2000) de- 
note the lines of flow, and that the spiral arms are trailing 
(this implies that the western part of the galaxy is closer 
to the observer than the eastern part). The resultant pa- 
ttern of projected velocities is shown in Fig.ll4H. Although 
the pattern is confined to the trajectory of the dust lanes, 
it agrees remarkably well with the projected residual ve- 
locity pattern in Figll4b in almost all respects. The only 
differences are seen towards the N and S of the projected 
inner disc, where the model predicts higher residual ve- 
locities than those shown in Fig. I14H . We explain this as 
follows: in the radial range from 10" to 20" our best rota- 



tion curve, shown in Fig. \5\ gives an apparent peak rising 
to some 300 kms -1 . It is probable that the peak, and the 
curve in this range, show enhanced apparent rotational 
velocity because we are not able to eliminate non-circular 
components in this range. Only the inner disc emits here, 
and opening up the free angle further reduces the signal 
to noise ratio to unacceptably low values. If we were able 
to make a rational interpolation between the rising rota- 
tion curve within 10" radius, and the slowly falling outer 
curve beyond 25", we might have attempted a more accu- 
rate separation of circular velocity and non-circular resi- 
dual velocity in the outer part of the inner disc. However 
such an interpolation would be so far from linearity that 
it would be little more than guesswork. Although the Hi 
distribution does not have the gaps of our Ha image, the 
angular resolution in the H I rotation curve given in Regan 
en al. (1996) is no better than 13" so cannot supply the in- 
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Fig. 14. a) iJST-NICMOS if band image of the central part of NGC 1530 with Ha isointensity contours overlaid 
at 1" resolution, b) Grey-scale representations of the J-K colour index image. Redder colours are indicated by the 
darker shades. Image taken from fig. 2 of Perez-Ramirez et al. (2000). c) Residual velocity map of 1" resolution of 
the central part of NGC 1530. This map was obtained after subtracting the model velocity field obtained from the 
modified rotation curve in the center from the velocity map (see text details). Superposed are isointensity contours 
of the unsharp mask from the iTST-NICMOS image shown in Fig. d) Projected velocity field in the line of sight 
assuming a gas flow of constant amplitude (60 kms -1 ) along the observed dust lanes (Fig. 114b) in the inner part of 
NGC 1530. All images, a to d, show the same area of the centre of NGC 1530. 



formation we need. The projected model field in Fig. I14H 
gives such good agreement over most of the spiral flow 
lines that we feel secure in claiming that the discrepancies 
to the north and south of the inner disc are most probably 
due to the effects of rotation curve subtraction described 
here. 



5.3. A scenario for circumnuclear spirals 

A generic explanation for the presence of circumnuclear 
spirals was given by Englmaier & Shlosman (2000). They 
showed that in a galaxy where there is a clearly identified 
ILR, the response of gas to forcing by a bar can give rise, 
for an appropriate central mass distribution and range of 
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Fig. 15. Unsharp masked NICMOS image 
(Fig. IHO of the central part of NGC 1530. 



sound speed, to a spiral pattern in the gas which can be 
present all the way down to within tens of parsec of the 
nucleus. They predicted that the arm/interarm contrast 
should be low, as found in near IR observations of these 
structures by Laine et al. (1999), and by Perez-Ramirez 
et al. (2000). The general framework for this kind of spiral 
structure does appear to be present in NGC 1530, since 
the form of the rotation curve has been shown by both 
Reynaud & Downes (1997) and by Regan et al. (1996) to 
give rise to a predicted ILR. Reynaud & Downes (1997) 
predict an outer ILR at ~ 1.2 kpc from the nucleus, and 
an inner ILR some 100 pc from the nucleus. Here all we 
wish to assert is the consistency of our results with this 
type of models. We do not attempt to extend the work by 
any type of detailed modeling exercise. 

6. Discussion and concluding remarks 

6.1. Three kinematic regimes 

NGC 1530 is a dramatically barred galaxy with very well 
defined morphological features, so it is not surprising that 
there is considerable recent work in the literature on its 
gas kinematics. However, the use of an Ha emission line 
map at ~1" resolution over the major part of the disc has 
enabled us to probe more deeply some of its kinematical 
properties. We have obtained a very well resolved rota- 
tion curve which is symmetric for the two independent 
halves of the galaxy. It was easy to detect the presence of 
non-circular components of the internal motion, in three 
dynamically distinct zones: streaming motions in the spi- 
ral arms, highly elliptical orbits around the bar, and an 
inwardly spiralling component particularly marked in the 
circumnuclear disc within 2 kpc of the centre. By the con- 
ventional technique of subtracting off a two-dimensional 
projected model rotation curve from the observed velocity 



field we produced a map of these non-circular residual ve- 
locities, which include components of amplitude well over 
100 kms -1 in the bar-induced elliptical orbits. The fact 
that the different kinematical regimes referred to above 
occur in distinct ranges of galactocentric radius, facilita- 
ted the separation of the circular and the radial compo- 
nents. However both in the range of the spiral arms, from 
63" to 90" from the centre, and in the inner disc, out 
to some 5" from the centre, we used auxiliary dynamical 
reasoning to effect an improved separation. We assumed 
that the departures from smooth radial dependence found 
in these ranges cannot be due to rapidly radially varying 
mass distributions, since no such variations are present in 
near-IR broad band images which trace the stellar mass, 
and must instead be due to projected components of lo- 
cal non-circular motions not well separated out from the 
true rotation curve. Based on this, we interpolated smooth 
underlying envelopes to the two sections of the observed 
rotation curve, and used this smoothed version of the ro- 
tation curve to produce revised maps of the residual ve- 
locity field, in which the streaming motions in the arms 
and the spiralling motions towards the centre in the cir- 
cumnuclear disc are better represented. These procedures 
do not affect the strong residual velocities around the bar, 
which occur in an intermediate range of radii, but do yield 
a more accurate overall map of the non-circular velocity 
field. 

6.2. The key role of velocity gradients 

By far the most notable new results in this paper refer to 
the gradients in velocity which we have been able to detect 
directly and measure. Particularly striking is the map of 
velocity gradients perpendicular to the bar, in the range 
of the bar (Fig. IlltO. Although in previous work on this 
and other barred galaxies the presence of a rapid jump 
in radial velocities, from ^+100 kms -1 to ~-100 kms -1 
as one transits the bar, has been detected, this is the first 
case where the relevant field has been closely mapped. It is 
essentially a shear field, since when we compare this with 
the velocity gradient map in the direction along the bar, 
(i.e. along the direction of the gas flow vectors) the per- 
pendicular gradient is considerably stronger. We do note 
that a perfect separation into components parallel to and 
perpendicular to the flow is not possible even in princi- 
ple in observations purely of radial velocities, but that the 
dynamics of the situation, with the x\ orbits around the 
outer lengths of the bar, is very favourable for a rather 
good practical separation in NGC 1530. It is notable that 
the zones of strongest gradient coincide precisely with the 
dust lanes even in fine detail. We have detected gradients 
as steep as 0.35 kms -1 pc -1 in a direction which is essen- 
tially perpendicular to the flow vector along the bar, i.e. 
these are shear gradients. Comparing this map with the 
Ha integrated intensity map in Fig. Illb we can see that 
there is a complete anticorrelation between high values of 
shear gradient and the local star formation rate. Massive 



18 



Zurita et al.: Ionized gas kinematics and massive star formation in NGC 1530 



stars are not being formed in zones of high shear. This 
is not a surprising result at all. The work of Rcynaud 
& Downes (1998) was directed exactly at this point in 
NGC 1530, and came to this conclusion. 

6.3. Parallel shocks inhibit star formation; 
perpendicular shocks stimulate it 

The only advance we claim here is to have completed 
and detailed maps of the phenomenology, with measured 
velocities and velocity gradients, both in amplitude and 
direction. From these we can extract several interesting 
semi-quantitative conclusions. Firstly, although star for- 
mation region avoid strong shear, they do not necessar- 
ily avoid high non-circular velocity. Specifically there is a 
star forming region, 20" to the NW of the galactic centre, 
in a region of high velocity along the bar (greater than 
100 kms -1 in a negative sense, deep blue in the residual 
velocity map). This region is, however, in a trough of ve- 
locity gradient: less than 0.05 kms -1 pc _1 both parallel 
to and perpendicular to the bar. It is interesting to note 
that this Hn region lies along the trailing edge of the bar, 
and that this edge is the site of a number of regions, while 
the leading edge of the bar does not show any strong re- 
gions. This is not a common situation and could lead to a 
conclusion which questions the placing of corotation near 
the end of the bar. However, we do not believe that this is 
necessary. The bar is so strong that the spray effect illus- 
trated in Fig. 03 should stimulate star formation along the 
trailing edge at positions of the arrow tips, where gas on 
quasi-elliptical orbits collides with gas in the bar as it be- 
gins to reverse its direction. Another way to consider this, 
is that we may be observing effects of a 4:1 ultraharmonic 
resonance, but highly flattened by the strong bar. All the 
other major star forming regions in the central disc, along 
the bar, at the ends of the bar, and in the arms, occur 
in regions of low residual velocity as well as low velocity 
gradient. 

Secondly, there does appear to be a relationship 
between shocks (as picked out by steep velocity gradients 
along the direction of the gas flow) and star forming zones. 
As seen particularly well in the map of velocity gradients 
parallel to the bar fFig. II lb) . linear zones of high veloc- 
ity gradient along the bar, at the edges of the central disc, 
and even in the spiral arms (not shown in the cited figure) 
tend to accompany the regions of strongest star formation, 
with the zone of steep gradient offset from the centre of the 
nearby star forming region by a distance of a few hundred 
parsecs. This phenomenon could have two causes: either 
the compression along the line of flow is causing major 
gas condensation which leads to star formation, or out- 
flow from star forming regions is impinging on the gen- 
eral flow to yield high gradients along the flow direction. 
We have considerable evidence from detailed kinematical 
studies of individual luminous Hn regions in a number 
of galaxies (Relaho et al. 2003a, 2003b) that outflow at 
velocities of between 50 kms -1 and 100 kms -1 is occur- 



ring in most of them, so the second phenomenon could 
well occur. However, the scales of some of the gradients 
seen in Fig. II lh and the fact that in some cases these oc- 
cur offset to one side of a star forming region and not to 
the other indicate that pre-stellar compression shocks are 
also being observed here. Quantification and modelling are 
clearly beyond the scope of the present paper, but these 
data may well offer a possible valuable route to the speci- 
fication of star formation conditions. It is, of course, the 
difficulties of pinning down observationally the conditions 
prior to star formation rather than any lack of theoretical 
offer which limits our understanding of the phenomenon. 

6.4. Quantified spiral inflow towards the nucleus 

A two-dimensional map of the non-circular velocities also 
gives us the chance to follow kinematically the gas flow 
which is picked out morphologically by the dust lanes. 
Although projection effects are quite severe on observa- 
tions of a flow whose vector in the plane of the galaxy 
is swinging through more than 270°, we have prima fa- 
cie evidence that we are able to follow this flow down to 
within 300 pc of the centre in this galaxy. This is not very 
close to the nucleus, but bearing in mind that NGC 1530 
is 37 Mpc from us we could certainly expect interesting 
results from this method for closer objects. 

Acknowledgements. We acknowledge Dr. M. Regan for kindly 
providing us his Ha rotation curve data, Dr. D. Perez-Ramirez 
for help with observations and Dr. S. Laine for directing our 
attention to the material in Figs, ll5l and ll4h and for discussions 
on its intepretation. We thank AURA for the use of the image 
of NGC 1530 in Fig. Et- The WHT is operated on the is- 
land of La Palma by the Isaac Newton Group in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto 
de Astrofi'sica de Canarias. Based on observations with the 
NASA/ESA Hubble Space Telescope, obtained from the data 
Archive at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. This 
work was supported by the Spanish DGES (Direccion General 
de Ensenanza Superior) via Grants PB91-0525, PB94-1107 and 
PB97-0219 and by the Ministry of Science and Technology via 
Grant AYA2001-0435. 

References 

Athanassoula, E. 1984, Phys. Rep. 114, 319 
Athanassoula, E. 1992, MNRAS 259, 345 
Begeman, K. 1989, A&A 223, 47 

Buta R, 1984, Ph.D. Thesis University of Texas, Austin 

Buta R, Block, D.L. & Knapen, J.H. 2003, AJ, (In press) 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H.G., Buta, 

R.J., Paturel, G., & Fouque, P. 1991, Third Reference 

Catalogue of Bright Galaxies (RC3) , Springer, New York 
Downes, D., Reynaud, D. Solomon, P.M., & Radford, S.J.E. 

1996, ApJ 461, 186 
Englmaier, P., & Shlosman, I. 2000, ApJ, 528, 677 
Elmegreen, B.G., Palous, J. & Ehlerova, S. 2002, MNRAS, 334, 

693 

Erwin, P., & Sparke, L.S. 2002, AJ, 124, 65 



Zurita et al.: Ionized gas kinematics and massive star formation in NGC 1530 



19 



Freedman, W.L., Madore, B.F., Gibson, B.K., Ferrarese, L., 
Kelson, D., Sakai, S., Mould, J.R., Kennicutt, R.C.Jr., 
Ford, H.C., Graham, J.A., Huchra, J.P., Hughes, S.M.G., 
Illingworth, G.D., Macri, L.M., & Stetson, P.B. 2001, ApJ, 
553, 47 

Greve, A., Reynaud, D., & Downes, D. 1999, A&A, 348, 394 
Huntley, J.M., Sanders, R.H., & Roberts, W.W.Jr. 1978, ApJ, 
221, 521 

Jenkins, A.A., & Binney, J. 1994, MNRAS, 270, 703 

Kennicutt, J.R. 1989, ApJ, 344, 685 

Knapen J.H. 1998, MNRAS 297, 255 

Lindblad, P.O., & Jorsater, S. 1981, A&A, 97, 56L 

Laine, S., Knapen, J. H., Perez-Ramirez, D., Doyon, R., & 

Nadeau, D. 1999, MNRAS, 302, L3 
Perez-Ramirez, D., Knapen, J.H., Peletier, R.F., Laine, S., 

Doyon, R., Nadeau, D. 2000, MNRAS, 317, 234 
Relano, M., Beckman, J.E., Rozas, M. & Zurita, A. 2003a, Rev. 

Mex. Astron. Astrofis., 15, 205 
Relano, M. et. al. 2003b, (in preparation) 
Reynaud, D., & Downes, D. 1997, A&A, 319, 737 
Reynaud, D., & Downes, D. 1998, A&A, 337, 671 
Reynaud, D., & Downes, D. 1999, A&A, 347, 37 
Regan, M., Teuben, P., Vogel, S., & van der Hulst, T. 1996, 

AJ, 112, 2549 

Regan, M., Vogel, S., & Teuben, P. 1995, ApJ, 449, 576 
Regan, M., Vogel, S., & Teuben, P. 1997, ApJ, 482, L143 
Roberts, W.W.Jr., Huntley, J.M., & van Albada, G.D. 1979, 
ApJ, 233, 67 

Rozas, M., Sabalisck, N., Beckman, J.E., & Knapen, J.H. 1998, 
A&A, 338, 15 

Rozas, M., Zurita, A., Beckman, J. E., & Perez, D. 2000, 

A&ASS, 142, 259 
Toomre, A., 1964, ApJ, 139, 121 

Tully, R.B., Fisher, J.R. 1988, Catalog of Nearby Galaxies, 

Cambridge University Press 
Shlosman, I., Frank, J., & Begelman, M.C. 1989, Nature, 338, 

45 

Young, J.S., Xie, S., Kenney, J.D.P., & Rice, W.L. 1989, ApJS, 
70, 699 

van der Kruit, P., & Shostak, G.S. 1982, A&A 105, 351 
Wray, J.D. 1988, The Color Atlas of Galaxies, Cambridge 
University Press 



